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ABSTRACT 


This report lays the groundwork for further development 
of a microcomputer that provides full celestial navigation 
capability. The physical design specifications and design 
philosophy are investigated. The fix algorithm which 
allows determination of position without the use of H.O. 
Publications is implemented in the BASIC language. 

The application of microcomputer technology to devices 
of this nature brings to light the need for an information 
distribution system. Thoughts on a system of this purpose 


are presented. 
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I. INTRODUCTION 


This paper deals with two interleaved subjects. The 
rwoO Subjects wilh be discussed with@varying levels of 
@etail and scope. The more specific subject of the two is 
the conceptual design of a celestial navigation computer. 

A general approach toward a digital information distribu- 
tion system is the second subject. The celestial navigation 
computer is discussed first, as it 1S an element of the 
entire system, and the design specifics of this device 
relate better to the entire system than vice versa. 

The celestial navigation computer provides for manual 
input of information as well as documentation of the out- 
put for legal records. The main reason why a computer is 
beneficial for celeStial positioning is that it eliminates 
errors made in recording, computing and plotting of celes- 
tial information. The design efforts attempt to eliminate 
manual recording of substantive information and manual 
plotting of lines of position (LOP's) for use in determining 
positdon. Automatic information recording is provided. 

This paper provides an analysis of the programming 
language, program I/O requirements and physiological require- 
ments of the device. The initial purpose of this paper was 
to design the device, then to see how this philosophy and 
design process could be applied to other building block 
devices applicable to the information distribution system. 


Comments to this point will be made in Section VI. 





Ti. NEED FOR A CELESTIAL MICROCOMPUTER PERIPHERAL 


It would seem that with all of the fancy new gadgetry 
on the market today, the age-old method of navigating by 
the stars might just start to fade away. It is conceivable 
Phat this could happen in the civilian world; however, the 
melrcary 1S Often more demanding inthe constraints placed 
Upem 1ts shaps. The electronic devices which use electro- 
magnetic radiation in the position determination loop are 
Subject tO deception, jaliming and detection (as in the 
ease Of trader). and support farlure, Geen , destruction of 
fiaviiedtronal satellites or EORAN stations, etc. The only 
other navigation device not subject to standoff disabling 
is the ship's inertial navigation system. It must be 
pointed out that an inertial navigation system is only as 
good as its platform's stability and the stability only 
meaecood as 105 position westimates. Thus the lomger the 
platform goes without being updated, the worse the platform 
memorlity becomes--truly a spiraling decline in accuracy. 

Pelle all Ofwehe electronic fiming devices ‘omt of the 
picture, viswal Sighting or sonar positioning is all that 
is left. Viswel sighting is passive and includes celestial 
DOsrelOn CeeemlirtroOn, aid ties 15 tite reason LOr Magn - 
t#ininge excellence in position determination by the stars. 

The mlerocomputer is a néwcomer to this field. This 


technology properly applied allows for increased speed and 





iecumae, and the reduction of human errors. <A need to 
determine position faster simply goes along with the increase 
mieonenemspeed, e€.g., mew hydrofoil ships can conceivably 
reach speeds in excess of 60 knots. The need for increased 
accuracy iS an economic as well as operational benefit. 
Accurate positioning is meeded for proper and efficient 
deployment of weapons, efficient navigation and maneuvering. 

The submarine force has implemented a celestial naviga- 
tion system in the most complete terms. Completely computer 
eontrolléd, seweral “spet shots"’ are obtained in a moment 
and the submarine's navigation systems can then be updated. 
The Air Force has used the ASTRO-TRACKER to aid in ceiestial 
Mavi? a tao. 

In January, 1975, Carl Nuese and Dr. Alan Schneider 
provided the algorithms necessary to create a celestial 
Navigatvoem™calculator usa@mg the Hewlett Packard Model 65 
programmable microcomputer [Ref. 8]. The peripheral device 
under design in this paper was intended to hit the mark 
somewhere between the HP-65 calculator and the submarine 
celestial navigation computer--much closer to the HP-65 
Side of the spectrum. 

The HP-65 calculator solution to navigation problem 
lacks the same sort of smoothness that placing a square 
block into a round hole does; with a lot of effort it works 
but it leaves a lot to be desired. The HP-65 general purpose 


calculator uses programmable magnetic cards to hold prestored 





programs that are executable on the machine. Though the 


device is relatively inexpensive the following drawbacks 


are noted: 


1. Storage and maintenance of the program cards are 
cumbersome. 

fm 2/0 -1S* nOnexis tent. 

5 Patri OrMa telson CONVeENtLOMNdl, 

4. Record keéping must still be done manually. 

ye) plne device 1S subject to. pilferage: 

Ge eEepanclon GCs Programs @aSedatiicult io ats 1d00- 


instruction memory. 


The advantages are the following: 


ie 


> WN 


Self contained (excluding the magnetic cards). 
Easier to Operate than most terminals. 
Inexpensive. 

Good Supportime Vendor. 


Applicable to many other mathematical problems. 
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Pirie DESIGN PHPEOSCPHY 


ch OESECTIVES 

The objectives of the design of the device are simple: 

1. Implement Nuese and Schneider's algorithms in a 

device that does not have the disadvantages of the HP-65. 

2. Maintain the advantages of the HP-65 without appre- 

elably Inereasine the cost/perftormance factor. 

See hovide fOr jimteeration into the overall information 

espn outl On sys tem. 

4. Provide the fail safe/soft features and hardcopy 

printout capabilities. 

Two definitions will be stated at this time. 
imenoceanpucecr: A small general purpose computer. i [t 
contains a microprocessor and a small number of (sup- 
port) LSI circuits. Usually mounted on one or two 


pr tnecd ee LiGiitanDOamaS  seiay OF May not be programnabile. 


Minicomputer: A computer with larger word sizes and 
faSter exeecitwen times than a adlerocomputer: 


The proposal to specify a microcomputer for the celestial 
Gomputer was based on the further observations that the mic¢ro- 
eomnce: can be exactly specified as to the I/O structure, 
Memom Structure amd contemts, power supply despens, ete, 
whereas the minicomputer requires adjusting of all of the 
above mentioned properties to "fit" the minicomputer. The 
cost for a microcomputer is also much iess than that of the 
minicomputer. 

An eight-bit-wide microcomputer consisting of a micro- 


prmecessor,ea l26 x 8 bit RAM@ea 1024 Sens ROM and peripheral 
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adapter cost around $473 in 1974, and by 1976 is estimated 
to cost under $100 [Ref. 6]. Though larger than the HP-65 
it is still less than ten times as large. An examination 
of data equipment suppliers wiil show that a small mini- 
computer will cost around $5,000 and higher. 

The device must be self-contained (except for power), 
easy to operate and relatively inexpensive. The I/O struc- 
ture must allow for keyboard entry, and the formats for data 
mucttepe: naturalized {(S60° 15" 33" vice 36.1553 ) to the prob- 
lem. Note that the algorithm expressed in Appendix A uses 
the latter form due to limitations of the BASIC language 
I/O structure. As the device would be a small production 
lot, vendor support would be minimal or nonexistent; yet 
the simplicity of the device should allow for easy onboard 


maintenance. 


Be  PHYSIGR DESIGN SPECI ELEGS 

The device must be housed in one unit. The unit must 
allow for connection to power and to the data distribution 
system. The unit will have a simple ten-digit keyboard with 
a N/E (north or east) and S/W (south or west) key. A 16 
element five-by-seven dot LED must be placed above the key- 
board for verification of entries as the line printer prints 
a COMM LCt camer tay, 

A Sgehnamaetem widestinesprinter (for printing instruc- 
POM moines ana cCaleilatron results} must be placed in 
the Uneeetneemnlt 1S 1lMistmated Invkueure 1, and a 


Sai eOWEptestray 1S displayed in Appendix B. 


ee 





Data Verification Readout (LED) a 


30 Character Line Printer 





NN . Program Selectors 


Data Entry Keyboard 


Fiseure 1-Celestial Navigation Microcomputer 


iS 





IV. THE CELESTIAL FIX ALGORITHM 


The purpose of the proposed navigational computer is 
to enhance the speed, accuracy and ease with which calcula- 
tions such as those used in celestial navigation can be made. 
The program presented in Appendix A removes the need for the 
user to depend on the Nautical Almanac [Ref. 13] and other 
meduinred H.O—m publications by computing the values supplied 
byeetumese tables. 

When taking a £1xX Several phases of sight reduction are 
considered: insertion of user values, correction of the 
Sextant altitude, establishment of the position of the 
observed celestial body, solution of the navigational 
triangle and displaying the results. 

User values include Greenwich date, day, month and year, 
the observer's assumed position, the name of the body ob- 
served, the time (Greenwich Mean Time), the sextant altitude, 
the observer's height of eye or dip and the sextant's index 
Bt Gorrection. 
| GCeomrection Of the sextant'’sS altitude iunceludes changes 
in altitude due to refraction, semi-diameter of the sun and 
moon and parallax for the moon. 

Establishment of the position of the obServed celestial 
body Onwtne cartnes(GP) amnvolves the calculatzon of the 
Greenwich Hour Angle (GHA) and declination (Dec) of the body 


(sec (re. 2) nscseOclatcd Wren Cac Stan ware eight 
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\ Greenwich leridan 


oy 
or 





Figure 2-Geographic Position (GP) or a Star 


as 





intermediate constants which are calculated and combined 
Vache Universal Constants {rate of precess#om ct the equinox, 
rotation of the earth on its axis and aberration) to produce 
the GHA and declination. 

Constants which apply to all stars are the following: 

ieee race of: rotation o£ whe earth with respect to the 
vernal equinox [Ref. 11]. 

2. The longitude rate of the sun based on the time 
between successive crossings by the sun of the vernal 
equanox [Ref. 10). 

5 wlihememisauity forsepoch 197S5maiRet. 10). 

4. The annual precession rates for epoch 1975 in right 
ascension and declination [Ref. 10]. 

With these common constants, particular intermediate 
constants associated with each star are developed. Utiliz- 
ing the apparent right ascension and declination as found 
in Apparent Places of the Fundamental Stars [Ref. 14], 
Fipht ascension and declination can be computed to a 
fetiemence time tineugh theweriectssofeprecession and 
aberration. Other constants include the rate of increase 
of the mean sun longitude in degrees per day, precession 
in right ascension and declination in degrees per day, true 
right ascensaoneand declination ame-degrees, and the time 
behavior of the corrections for aberration [Ref. 8}. 

When the intermediate constants are developed and 
applied to the algorithms, the computed GHA and declination 


of a star will result. Discussion of the computation of 
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eiewconstants, Ginvand declination iS presented in the 
next section. With the algorithm for computing GHA and 
declination of a star implemented it now becomes a matter 
of selecting those stars necessary to perform a Sight 
reduction yielding a navigational fix. 

Consigerine a thwee=smaan fix, the GHA and declination 
for each star is obtained as discussed above. Working with 
pairs of stars the solution of the associated navigational 
triangles will produce intersection points on the circles 
of equal altitudes (see Figure 3). 

Three triangles are used to calculate a possible position 
as shown in Figure 4, The first triangle has vertices N 
imontn pole), A and B. The angle A 1S either Gl1-G2 or 
360°-(G1-G2). This angle represents the difference in the 
GHA's. Given the angle A and cosides dl and d2 (the 
declinations), the law of cosines 1s applied to solve the 
fmind Side t. Given the three cosades the angle 8 opposite 
coside dl is calculated. Given three cosides in the second 
tueamele, Hi, HZ and t, the angle BS opposite H1 is found. 

Pceulne sie Deetie Sitio Oreelc San(Gl-G2Z), the angles 
B(1)=mB-B3 and B(2)=mB+B3 are found. Given angle B(1) and 
eosides a@2.and HZ, the third coside L{1) 1s found. Given 
the three cosides d2, H2 and L(1) the angle A(1) opposite 
221s eottuval, In the same manner L(2) and A(2) are found. 
L(i) is equivalent to the latitude and A(i) iS equivalent 
Commence lomerelce.O; Aaneintersection of the circles of equal 


altitude. 


Le 





Star A (GP) 


Circles of Sees 
Equal 90-H1 —s 
Altitude 7 


mE: 


90--H2 


Sta 8 (08) 


Figure 3-Navigational Fix from Altitudes 
of Two Celestial Rodies 
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N (north) 





Figure 4-Triangles used for Navigaticnal Fix 


bes, 





Repeating the algorithm considering the third celestial 
body and one of the previously used bodies, a second pair 
of positions will be found. A third repetition yields two 
more positions totaling six (the maximum) possible positions, 
as shown in Figure $§. 

These six positions are then investigated to determine 
the closest three by determining the shortest distance 
connecting these three points. With this resulting triangle 
moe location of the two endpoints of a line which is one- 
mitd aWway from any side of the triample is computed. The 
center of this line containing the two endpoints is the 
center of mass of the triangle and therefore is considered 
me celestial f1x. 

The program will provide for accuracy in mathematical 
manipulation of numbers to one part in four billion using 
mee DbASIC interpreter at the Naval Postgraduate School and 
miplementing the program on the INTEL 8080. This will pro- 
vide for the accuracies of .1 or .2 miles as presented in 


Ref. 8. 
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Figure 5-Overall Celestial Situation 
with six intersection points depicted 
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V. PROGRAM THEORY AND IMPLEMENTATION 


A. THEORY 

When computing the constants needed for the program 
a common reference time 1S necessary to provide for the 
@eccuracy ot results. if the data is referenced to the mean 
equinox and equator of date, or to a fixed epoch then 
no error is introduced. The reference time employed in 
the algorithms is 0000 29 Nov 1971 as suggested by Ref. 8. 
Any data known at some other time is retarded to this 
reference time. 

Note that the GHA of the star is related to the 
Greenwich Hour Angle of Aries (GHA) and the sidereal hour 
angle of the star (SHA), the latter being 360° minus the 


right ascension a. Thus for any time t, 


GHA = GHA, + SHA (1) 
GHA = GHA, + (360 - a) (Zz 
GHA = GHA, - oO ) 


The GHA, increases at a constant rate, namely, the rate 
of rotation, a, of the earth with respect to the vernal 


equinox. Knowing GHA, at some time ty, Giten 


eee a (cme GHA, | - a(t) (4) 


Lil 





where t is the present time measured in days and decimal 
parts of a day from the time ty, a time in the recent past. 


Because it will later prove efficient to refer all events 


to an earlier time Ger 0000 29 Nov 1971) this leads to 


GHA = a(t - tad - a(t - t,) ap GHA, 4 - a(t) (5) 
= a(t - eed F GHA, - a(t) (6) 
where EO = GHA - a(t, = to) (7) 


The apparent right ascension ao at time t is related to 


its known true value a' at some previous time t, through 


1 
the effect of precession and aberration. 


Therefore (see "Supplement," [Ref. 11, pp. 39-49]) 


a) To Gap c(t, ) Ao (z) (8) 
where 
a'(t,) = the true right ascension at time ty 
OW = precession in right ascension (degrees 
per day) 
Aa = COLRVeGmIOn In Tron aScenSion,at time mt 


for aberration (as defined on p. 4/7, 
"Supplement" ) 
In order to refer all data to the common reference time 


t it 1S convenient to rewrite (8) as follows: 


QO? 
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GE) = eer) a(t - a) a(t, oe enere ont) (9) 


ke mec CC ACNE) (10) 


where 
con e= Ge) a(t, oe) (mk) 


) 


imsecrtineg (10) into (6), 


GHA = GHA, "aCe = tJ - a'(t.) - a(t - tJ - hatte) ei) 
GH re Meee es Ca Ce te) sew Aa (tr) (13) 
Or, 
GHA a, * Sag tt — eae Aaitt) (14) 
where 
ass = GHA, Se UIGaR) (15) 
= CHAS - a(t, = a, 
Secu Geap ea (0G) Pate) (16) 
and 
ajg = a - a ey) 


Now Ao must be developed into a usable form. The equation 
On pages 47 and 48 in "Supplement" [Ref. 8] shows that aber- 


ration at time t is given by 


= Khpscewommicosua (cos Awcos « + Sin @’ sin A] (18) 


Where A is the longitude of the mean sun meaSured along the 
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ecliptic from y, e€ is the obliquity (angle between the mean 
planes of the ecliptic and equator), k is the constant of 
aberration defined on page 48 in "Supplement," having a 
mame Of 20.47", and ao" and 6' are the true right ascension 
and declination, all values being at time t. Selecting ty 
as a convenient time for determining the mean place of the 
Sear and the obliquity, and regarding a', 6' and € as 


@enstancs over the interval from t, to t, equation (18) 











1 
becomes: 
= AO OE 
Aa = - (sEG0 SeC 64 cos ory cos e,) COSex 
- (Sen SEC 64 Sin a7) Sin i’ (19) 
= a1 900s A + a4 Sin r (20 
where 
=" 20a 7 
a5 Z6007 Sec 64 COs Oy COS €4 PZ) 
2 Ona ae ; 
a54 ZEOH SES 64 Sin a, (22) 


The factor 3600 enters because Aa will be computed as 
fractions of a degree whereas k was given in arc seconds. 


Equation (20) can be rewritten as 


Aa = a Seaiien(cue ta (250) 


12 13) 
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a fae 


& -—r 


a 





where 


7 2 2 
a) bl Sa ea (24) 
and 
_ elite 150 
eae a eat 7? (25) 


The longitude of the mean sun increases linearly with 


time according to the equation 


NS Aen. t,) +A (26) 


it 


where dy 1s the known longitude at some convenient reference 


time t,- As has previously been done, it is best to refer 


to a common reference time t, by rewriting (26) as follows: 


Ga Con nc 


A(t) ak eae CAR 


hate ce) at Ne (28) 


where ro Vomdebined as follows: 


ae) (25)5 


Pepstitutine (28) into (23), the results are 


Aa =a C20) 


12 Sin [A(t - t,) - A + a 


i3! 


ay, sin [A(t - t.) + aqgl (31) 
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: = a - 
oa 7 

. = a y 7 7 - . : a ; we S a 
a _ ee A ee ee a : a _ 











where as 1s defined as follows: 


as4 7 Le <r a43 (iS2) 


Using (31) in (14) the substitution yields 


GHA = aig * ajyg(t - ty) - a5 sin Petey ip él (33) 


14! 


Because aj9 Inciudesmecdwen rateyeroughly 360° per day, and 
enat t - t, may be on the order of 30,000 days, to preserve 


accuracy in the computation it is worthwhile to separate 


249 into two parts, one of which, Ao» ts 5900 Wand the 
other aoo> 1s the remainder. Let 

1S peo 7 (54) 
where 

a5, = 360 (35) 
and 

a5 5 = 249 = Sow (36) 


In similar fashion, denote the number of days from the 
reference ty to the time t to be D.d where D is a whole 
number of days since 0000 Nov 29, 1971, and da fractional 


mart. That 1S 
fears (Co) 
Then it is obvious that the product of the angular rate times 


the time yields an angle modulo 360° which can be expressed 


oe 





aise o | lows: 


ayg(t - Ee) SS (et0) st aoo) (D cn) (38) 
= 360D + A549 (D dete oO Ud (39) 
= a4 (D jm ec 1 5600 (40) 
= a4. (D.d) + 360d (41) 


PupSsStituting this result and (37) into (33) yields 


GHA = + a> (D.d) 0 0 Gane a, 7[A(D.4d) + ay 4) (42) 


ai be 


The constants o and §, precession in right ascension and 
G@eciination, for a particular star are oDtained using the 
formulas on page 39 in "Supplement," except that the values 
of the constants have been taken for epoch 1975 from the 
American Ephemeris and Nautical Almanac [Ref. 10]. Conver- 


St1on to degrees per day 1s 


P 1 
® = Senn x 365.24719418 [46.1060 + 20.0404 sin ay tan 64] (43) 


Ove 
! 


c= Besa [20.0404 cos a, ] (44) 


This now completes the derivation for the Greenwich 
MOUr Angle: Or a Star at time t. 

Turning attention to the derivation of declination, it 
is noted that the derivation is similar to those just 
Garried out. Declination of a star changes because of 


prececSsmemmeand parallax. Letting 6 be the apparent 
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declination at time t, 6' be the true declination at some 
Specitied time, 6 be the rate of precession of this star 
in declination, AS bemthe Gorrection to declination for 


aberration, the result 15s 
6(t) = é'(t,) + 6(t - t,) Teno tees) Clay 


where ty 1s the reference time of which the true right 
ascension and declination are known. The numerical value 
of 6 is given in (44). This form is converted to the 


standard reference time t. by the same method which has 


been used above. The result 1s 
6(t) = ote) spe (6) (Oem t) rao ( t) (46) 


where 


g'(t) = 8'(t,) - é(t, - t) (47) 


From pages 47 and 48 in "Supplement" the time behavior of 


mae Correction for aberration is 


Moe — -k(Sin Gyncos 6!.- cos €, sin a! Sin $3) COiSas 


! 
it 





1 1 1 
-k cos or Sin 6F Sin i (48) 
= a). COs A + 246 Sin A (49) 
where 
a = 0.5! if(come. COSMO. .—- €OS O. Sain Oo. “Sin 6.) (50) 
ce 3600 es > | ! i 1 
_ 20.47 ae 
a5 6 ZEOH COS 4, Sin Sy C5?) 


Zo 





ine conovents a and ay¢ are i” degmees. clmplifyine further 


5 


a 
err i@ lector a) (52) 


Aé =a 
ley 24 6 


where 





2 Z 


1 Ste Ts (55) 


pepsi tutine (52) and (28) into (46), 


6Ge) = omGes) oad) + ai7 S leant ee t,) - aq! (54) 
where 
-1 “15 
a = ) eeetan — (55) 
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The final equations.of importance which lead to GHA and 
declination are (42) and (54). Note that eight star- 
particular constants are involved in these two equations, 
plus ie a constant which appears in all calculations of 
GHA and declination of the stars. 

The equations involving reference times have been set 
ip in such a way as to maximize freedom in their selection. 
Reference time to was selected to be 0000 Nov 29, 1971, to 
facilitate the calculation of time in days with a minimum 


number of entries as suggested by Ref. 8. 
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Messeteom@ence time ty 1S any convenient time, presumably 
Pe ehes recent past, atewhich the GHA of Aries is precisely 
Piven in some reference tabulation. The time ty is referenced 
to true right ascension and declination as listed in some 
reference almanac. The true right ascension and declination 
at a'(t,) and 6'(t,) can be obtained from the apparent 
values a(t,) and 6(t,) by setting t = ty in equations (8) 
ome, (45) and solving for c* and 6'. 

The reference time ty is also the time, again in the 
recent past, at which obliquity Ey: Oy > 54 and dy are 
accurately known. The reference time for hy is such that 
the accuracy of the sun's mean longitude is known from some 
oimanac data. It will minimize the errors due to (a) ecden- 
meeeity Of the carth's orbit, amd (by) the leap year cycle, 
if ty 1s choSen at or near the time of earth perihelion, 
mre tne Midpoint of a leap year cycle. 

It should be noted that ty for the known value of GHA 
of Aries may be different from the ty for the known values 
of true right ascension and declination. Likewise t, may 
be yet another time for the known sun's mean longitude. 


Since all the values (constants) associated with t, are 


retarded. to t, no ditterence in ty TS Ob jeetnonab le , 


B. IMPLEMENTATION 
The program presented in Appendix A is the result of 
several conversions. The first conversion consisted of 


taking algorithms expressed in Ref. 8 with the HP-65 
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programs to produce a debugged FORTRAN IV program. The 

HP-65 programs deal with degrees as a unit of angular measure 
vice radians. A conversion of all data from degrees to 
radians was made for use in the FORTRAN IV and BASIC pro- 
grams. General routines 0010 and 0020 provide for I/O 
conversion of radians to degrees, and hours, minutes, and 
seconds to hours (and vice versa). These conversions were 
needed for ease of data entry and data output. 

The HP-65 programs are compact due to the limited 
number of registers used for intermediate storage. This 
efficiency of memory was carried through to the FORTRAN IV 
programs and then to the BASIC programs for the routines 
that determine the GHA and declination of the star and the 
subroutine that determines the azimuth (Zn) and distance 
(subroutine 0300 and 0400, respectively). It was felt that 
this level of efficiency need not be continued in the 
remainder of the program due to the need for self-documentation. 

FORTRAN IV was selected as the intermediate step in the 
conversion due to the amount of core memory allocated to 
the FORTRAN compiler. The BASIC interpreter is limited in 
memory on the IBM 360 installation and would not hold the 
completed BASIC program. 

The reduced precision of the FORTRAN IV version of the 
program (as it was not programmed in double precision) was 
attributable to the internal formatting of the variables. 
High precision was not of great concern at this stage of 


development. 
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As the goal of this process was to provide a program 
that would implement the algorithms necessary to complete 
a celestial fix on a microprocessor, the FORTRAN IV version 
was converted to BASIC. The BASIC program allows for much 
greater precision in computation as the values are expressed 
in 32 bit floating point, when interpreted by the available 
meokC anterpreter in the Microcomputer Laboratory at the 
Naval Postgraduate School. The BASIC interpreter at the 
Wwe C. Church Center was not of sufficient size to run the 
entire program. Therefore, the main Subroutines and all 
required support routines were run and debugged separately, 
and required information was transferred manually from rou- 
tine to routine. Because the BASIC language does not pro- 
vide for scopal variables and due to the limited types of 
identifiers available much of the self-documentation of the 
FORTRAN IV program was lost. Table I correlates variables 
in subroutine 0100 in the program listing (Appendix A) to 
those variables previously discussed during the theory 
presentation. 

Subroutine 0100 computes the eight constants associated 
ween a Star and places them in array C (elements C(2)-C(9)). 
Element C(1) is the rate of increase of the mean sun longi- 
pide (commom to all stars). 

In order to organize the variables as much as possible 
a scheme was invoked as explained below. Since there were 
a number of subroutines in the program, a method of passing 


Viltaples was. conseructeda. The variables’ RI-R9 and P1I-P9S 
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TAB Eee! 
idemeutrer ConrmeLadwmon of Varaables ini Subrowtine 


0100 and Equation Variables 


Subroutine Variable Equation Variable 

C(1) A 

CHEZ) Ao 6 
GES) a47 
(4) 3 

C5) §'(t) 
Cina) a4 4 
C(7) as 
C(8) a1 9 
C(9) ay4 

is 04 

[4 $4 

15 3 

[6 as9 

{7 as9 

[9 a4 
JO arc 
Jl 246 
GS Ey 
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are passing variables. Variables used in three or more 
non-general routines were considered as global and were 
selected to be as self-documenting as possible. For example, 
the two constants for each star are placed in the array 
S(X,X). The constants used by routine 0300 to determine 

GHA and declination which are created from the S array are 
placed in the C vector. All other constants were repre- 
Sented by the format C1-C9. 

Several routines use intermediate variables while calcu- 
lating the object variables. These variables are repre- 
Sented Dy Li-19, J1-J9 and MI=M9. All logical values are 
represented by L1-L9. Note that this scheme was invoked 
to eliminate problems created by not having scopal varia- 
bles. This must be considered in the selection of a 
language for other peripherals used in the system. Points 
to be considered in selecting a final programming language 
for such a system and the peripherals attached to it are 


eascussed in Section VII. 
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Vie THE -INPORMATION DISTRIBUTION SYSTEM 


im  OpnvseGrives 

This subject will primarily be concerned with the 
integration of a general purpose computer into the existing 
shipboard information flow. The "system'' in the spotlight 
Momence iaviegation Suite aboard U.S. Navy shaps. The phil- 
osophy for this approach is simple and is best stated by 
Joseph Weir [Ref. 1]: 

"Good systems can be designed only in the presence of 
potential users and applications. By using the system 
in numerous applications and assiduously watching its 
performance, it gradually can be adopted to the (actual) 
needs of the application and the user. Ingenious ivory 
tower solutions frequently do not work or they solve 
pPeoolens that don't. exist)’ 

A close examination of this statement while relating 
it to the navigation suites will show that in fact the 
"system" has been designed both as far as it need and can 
be with the existing equipment. Note that the "system" 
comprises geographically dispersed equipment and telephones 
for communication links. 

The application and needs of the navigation system have 
been defined as will be discussed later in this section. If 
all this is true then why rock the boat by imposing a 
system, apparently complex, on an existing "system''? 

Discussion with crewmen aboard the USS Hancock indi- 


cated that the existing system is cumbersome, subject to 


human frailties and is not timely enough for the devices 
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on board. OMEGA radio navigation sets are capable of 
providing information for navigational fixes every few 
Scecomess, yet this information is only used hourly. In 
this case the communication link is the "weak" link, and 
to eliminate this weak link an information distribution 
System must be imposed. The distribution system must 
appear translucent and straightforward to the user, for as 
Weir mentions, 

bine more complex the system, the more likely it is 

that it won't get used or will be misunderstood or 

mesused.'' [Ref. 1] 

All systems have sensors and effectors, and some 
communication constraints placed on them. Although most 
of these sensors for the navigation system in fact exist, 
there are always new devices being created. 

For example, ten years ago operational OMEGA navigation 
Seesedid not exist. If a system had been created and imple- 
mented on the ships at that time without forethought with 
regard to new devices, then the information system would 
not have been able to integrate the new device into the 
network. 

ii Onde etomeuiy utilize aneineormation distribution 
System there must be "terminals" with which to input informa- 
tion and query the system. For example, suppose the bridge 
wWantS tO know the average RPM of the four propulsion shafts. 
Engineering would have a terminal (either automated or manual) 
GidtmerouldudececepeusGcertain required data, €.g., tachometer 


MipiecmebaistOrrm tC 1nto 1tS varied forms and relay the 
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data to the central processor. The processor would then 
transfer the required information to the terminal at the 


bridge. 


ie DESIGN PHILOSOPHY 

The question might be asked, "Why try to integrate an 
miormation distribution system with the existing periph- 
erals?'" The answer should be obvious--cost of implementa- 
tion. Existing equipment must be fully exploited until 
it 1S not wise to continue using it. 

By integrating existing equipment into a network instead 
of creating all of the component parts of the system anew, 
the cost of research and development of all the new sensors 
and effectors in the system can be elided [Ref. 2]. Of 
course, the use of current equipment will demand the design 
of interfacing equipment such as analog to digital converters, 
digital to analog converters, fault sensing networks, etc. 
However, the cost of these units should be relatively inex- 
pensive when compared to the redesign of any one sensor or 
effector in the system. The application of the same item 
(A/D or D/A, etc.) to many peripherals is a cost savings 
when design efforts are considered. 

ASWdiStrIbution of existing data is the concern of the 
central system and since there is a multiplicity of like 
data, the question of which source's information should be 
distributed must be handled. A rather trite answer is 
"the best."’ Again a question, "Which sensor is providing 


pene woes = anLoclatirvon: 
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Sade einitenen Of vthe term "best™ 15 recursive and must 
micwude WOrds lake timeliness, accuracy and statistical 
Valwadiy, ASSeden Sensor in the system Gs operating at 
pre-system design rates, some periodic (auto track LORAN C) 
and others at non-periodic rates (celestial sightings), 
timeliness takes on a meaning associated with the quality 
of the source. For example, if the only peripheral capable 
of input to the system is the celestial navigation unit, 
an output of ship's position would only be considered timely 
if the last celestial position had been updated with inputs 
from the "best'’ heading and speed information to provide an 
estimate of the current position. 

The word "accuracy" when associated with "best" means 
most accurate on a timely basis. A position based on a 
ten-hour old OMEGA position is usually not as accurate as 
a five-minute old celestial position. 

This 1S all well and good, yet how are these descriptors 
compiled? A central machine is obviously the device that 
will have to provide the software to meet problems based 
on statistical properties. 

Statistical properties of each measuring device must 
be taken into account. The proper error probabilities must 
be applied to each input; then based on the result of the 
combined statistics, the "best" source is selected to repre- 
Seme that classsef input. It is this source's information 


(corrected for time) that must be distributed. 
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Along with statistical constraints imposed on the data, 
a status check must also be imposed on the device. If 
Pic ddte@wis apparently correct yet the device has failed, 
then the data is obviously of no use. The status determina- 
tion system must be as nearly fail-safe as possible. Since 
mis Paper 1S proposing a Centralized star network vice a 
jmery Connected or ring network [Ref. 2], all the status 
checks must be initiated at the central computer (see 
Figure 6). In order to design a fail-safe status monitor, 
the following considerations must be employed: 

1. Notice of failure of the central computer must be 
arstributed to all peripheral stations. 

2. Tamlure Of each peripmenalemust affect only that 
peripheral and, at worst, place the system in a degraded 
mode. 

3. Status of selected peripherals might have to be 
recorded for legal records. 

4. Failure of the peripheral must not cause failure 
Of the status check mechanism. 

There are several reasons for designing an information 
distribution system. Many operations aboard a ship are 
dependent on the same source of information. For example, 
the Combat Information Center and radar stations need true 
course and speed information. These are parallel operations 
and their information must be based on the same source. A 


distribution system would be the answer to this problem. 
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All of the discussion thus far has been about distribution 
Gf information gathered from sensors, so where is all this 
fine information used? An example may best describe the 
answer. Suppose the navigator updates the current ship's 
position with the celestial navigation computer. If the 
Captain desires to proceed to a desired point, the point 
could be placed into the system and the system would pro- 
vide proper steering information (great circle steering) 
to the Helmsman, or this data could be transmitted to an 


aaeo-pilot upon request. 
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VII. LANGUAGE SELECTION 


The subject of language selection is placed at this 
point in the paper due to its applicability to both the 
celestial navigation microcomputer and an information 
distribution system as described within. The considerations 
in selecting a programming language for these two devices 
have much in common, the most important of which is the 
Microprocessor . 

There are two main considerations in selecting a pro- 
gramming language for a microcomputer. The first and most 
mimortant Consideration is that of efficiency; efficiency of 
both memory utilization and program execution speed. The 
second point concerns the cosmetics of the language (i.e., 
is the language easy to read and is it self-documenting?). 

imeoncideration of the first propere, we thateon 
efficiency, the language must be chosen in the light of 
cost/availability of volatile memory. Even though the 
addressing schemes of some microprocessors provide for 
large memories there may be physical and/or cost limita- 
tions placed on the microcomputer design that make very 
efficient use of the available memory imperative. For 
peripherals operating with unmodifiable code the number of 
memory locations can be optimized and placed with a high 
degree of efficiency. For processors using programs which 


have subroutines with arrays of data created during the 
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excettmion of the program, the placing Gi these arrays 
On aemetal mask chip or fixing and reserving the addresses 
of these arrays is an inefficient use of memory. 

The best way to program in an efficient manner is to 
use block structured programming techniques utilizing 
"scopal" variables. Block structured languages like Algol 
Seeand PLM (Ref. 12) allow for scopal variables to include 
arrays which provide for efficient use of random access 
memory and reduce the amount of programmable/read only 
memory. By creating arrays as needed, common random 
access memory may be used very efficiently allowing for 
lower investment in space and cost. 

The second area of concern, that of cosmetics, 1S a 
growing concern with all languages. It may be a seemingly 
unimportant aspect until the prospects of non-volatile 
runtime modifiable LSI memories (e.g., ferroelectric 
memories [Ref. 15]), multiplicity of similarly constructed 
peripherals and common maintenance support are considered. 
The celestial navigation microcomputer program bases 
ait@f the calculations for geographical positioning 
of the body constants which vary from year to year. 

In order for the device to enjoy extended longevity, 

the constants must be alterable. The language must 
allow for the documentation to help provide this support. 
An assembler program is not a self-documenting 


program no matter how efficient it is. Its logic and 
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meaning are often forgotten with the loss of the programmer. 
Indeed, the programmer may lose the logic to his own program 
if he has been away from it for a length of time. For 


example, the statement below is obvious in its intent: 
NO$YEARS = NO$DAYS / 365. 


Yet even though the next statement provides for the same 
Smerelency Of calcnlation it does not “eonnote as clear a 


meaning: 
Pio. Gl. 5605: Dive LIT ye: solO . 35. 


The assembler language form of this statement loses the 
easily conceived meaning of the first statement. 

Though the celestial fix program in Appendix A was 
written in the BASIC language due to the precision attainable 
through its use, the program itself would not be under- 
standable without the comment statements and extra 
documentation. 

Specifically, there are primarily two high level 
languages available to the microprocessor programmer ; 

BASIC as implemented at the Naval Postgraduate School and 
PLM. BASIC, as the name implies, provides for the basic 
needs of a programmer. BASIC has subroutines, iterative 

and logical capability, the full range of mathematical 
manipulators and limited I/O. A quick glance at the program 
in Appendix A provides the evidence of a major drawback: 


Ldeke ome varbtable scope and identifiers. — There ane moulorad 
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variables; all variables are global. This property coupled 
with the limited identifier schemes causes a programmer of 
ieee program to run Out sof variables quickly. His atten- 
tion is turned to efficient use of identifiers instead 
of efficient programming techniques. 

PEM) On the other hand, provides for sceopal variables 
and identifiers that allow for self-documentation. The 
following statements are equivalent; however, the second 


connotes more meaning: 
LET Ji = G(Q) - G(R) 
DIFF = GHAl - GHA2. 


The array used in the determination of the three 
intersections used to calculate the fix (subroutine 0600 
ARRAY V(6,6)) 1S a temporary array. In BASIC the erray is 
permanently allocated space. Using PLM and scopal variables 
the array would be created by the program during execution 
(indeed, cculd be of variable size) and be destroyed or 
written over upon exit from the subroutine. There are a 
few drawbacks with the use of PLM. The primary drawvack 
in the execution of the program in Appendix A is the wora 
Size. Sixteen bits do not provide the accuracy needed for 
me celestial navigation problem. Though this problem can 
be solved, it is not easily solved. 

Special mathematical routines like trigonometric func- 
fLONS Are nNoteavallable in PLM, yet they can be programnica 


iiempeune 1S that PLM can be tailored to tne meéeds on Femme 
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pe anor With a little extra work while the BASIC language 
Cannot be tailored to the user at ail. PLM provides a 

much more versatile language yet the detail required in 
programming is greater than that of BASIC. 

A final comment on the selection of a language deals 
with the multiplicity of processors on an information 
distribution system. If each processor were programmed at 
a different level with a different language the maintenance 
and design problems would be horrible. For example, suppose 
one processor is to provide a second processor information. 
If the first processor is programmed in BASIC then it would 
probably transmit its data in 32 bit floating point form. 
The second processor (programmed in PLM) would have to have 
an I/O routine to convert that processed data to a local 
format. If there were a third language then the PLM- 
programmed processor may need yet another [/0O processing 
routine. This Scouaclon 1s the grounds stor se Leectlonsorea 


single programming language for a system of microprocessors. 
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Val. SICONCLUSTIONS 


Currently, the proposed computer is capable of developing 
a navigational fix (without use of the Nautical Almanac and 
feo. Fublications) utilizing only the 57 navigaeronal Stars. 
Because of this limitation it is obvious that numerous and 
known applications are yet to be incorporated. Among these 
include methods which employ the attributes of the sun, moon 
and planets so that positioning with these bodies can be 
performed. Considering functional aspects for future 
development of programs, suggested programs are closest 
Hount Of approach (CPA), great circle calculations, star 
location utilizing estimated altitude and azimuth and £f1x 
er langulation routines. 

The sun, moon and planets present special problems 
because of their proximity to earth. They do not appear 
fixed in the sky as do stars. Stars do move across the 
Sky due to the earth's rotation but their great distances 
impart a constant relationship to their neighboring stars. 
The planets and moon, being closer to earth, possess no 
fixed relation in the evening sky with the stars. 

The size of the sun and moon requires a special con- 
Sideration termed semi-diameter. More accurate sight reduc- 
tions can be performed by aiming the sextant at the upper 
and lower limb of these bodies; corrections for the radius 


of the body can be programmed. 
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For lunar observations horizontal parallax is a 
parameter required to be considered. It is the apparent 
difference in the position of a body when viewed from two 
@itterent locations. Horizontal parallax is significant 
for the moon because it is the closest body to earth. 

The celestial lines of position are included to give 
me NaViegavor the option of plotting LOP's. This method 
Should provide the navigator with an azimuth, a distance 
and a toward/away indication in reference to the assumed 
position of the observer with which the LOP should be 
petted. Thais 1s the backup routine to the {1x calculation 
algorithm. 

The proposed celestial navigation computer was used as 
a vehicle for the approach to a design of a dedicated micro- 
computer. Other devices of similar construction used as 
peripherals in the information distribution System would be 
approached in the same manner. Approaching the design of 
each device on its own merits will allow for the system to 
pEow in a manner that "fits" the situation, i eC ces ime le 
mivory tower’ system would not result. 

Indeed, the development of the central processor 
(assuming enough ports are planned for) would grow with the 
system as each new unit would simply add length to the polling 
routines and interrupt structures. Following along with 
Weir's natural growth concept, each processor's programs 
would be designed and built (or at least the algorithms be 


provided for) as needed by cognizant individuals, from 





Glierene Laterature, as done in the case of the celestial 
Maverweation computer. The only restraint placed on the 
design/designer would be that of language and interface 
ceri tedtions, As delineated carlier, the authors believe 
that of the two languages examined PLM would be the most 
efficient language and would provide for much needed 
pele Gocumentat lon. 

ine seceneraia processor Ian the information <tst1 1 buco 
System may very well be a minicomputer or a microcomputer 
fepending On 1ts functions. it need only be a distribution 
device designed along the same lines as the peripherals 
mee supports. the level of computation provided by the 
central processor must be established early in the design 
phase because the peripherals would need to "take up the 
Slack" for the computational capabilities not provided by 
the central processor. For example, the averaging of the 
revolutions per minute may be accomplished at the peripheral 
(which is recommended) or at the central processor on the 
request of the inquirer. Again, keeping this device 
Simple (only the necessary switching and monitoring routines) 
will provide for more independence and faster information 


transter. 
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APPENDIX A 
GERESTIAL FPUX PROGRAMELISTING 
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